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Abstract. By extending the breadth first search algorithm to any d-type critical or 
subcritical irreducible branching tree, we show that such trees may be encoded through d 
independent, integer valued, d-dimensional random walks. An application of this coding 
together with a multivariate extension of the Ballot Theorem allow us to give an explicit 
form of the law of the total progeny, jointly with the number of subtrees of each type, 
in terms of the offspring distribution of the branching process. 



1. Introduction 

Let Ui,U2 ■ ■ ■ be the enumeration in the breadth first search order of the vertices of a 
critical or subcritical branching forest with progeny distribution u. Call p{ui), the size of 
the progeny of the i-th vertex, then the stochastic process {Xn)n>o defined by, 

Xo = and X^+i - = p{un+i) - 1 , n > 

is a downward skip free random walk with step distribution F{Xi = n) = v{n + 1), from 
which the entire structure of the original branching forest can be recovered. We will refer 
to this random walk as the Lukasiewicz- Harris coding path of the branching forest, see 
Section 6 of [Tj, Section 1.1 of [3] or Section 6.2 of |13) . A nice example of application 
of this coding is that the total progeny of the k first trees ti, t2, . . . , tjt of the forest, see 
Figure [H may be expressed as the first passage time of (X„)„>o at level —k, that is, 

Tk = inf{z : Xi = -k} . 

This result combined with the following Kemperman's identity (also known as the Ballot 
Theorem, see Lemma 5 in [21 or Section 6.2 in 1131): 



F(T, = n) = = -k) , 

n 

allows us to compute the law of the total progeny of ti, t2, . . . , in terms of the progeny 
distribution v. Note that the total progeny is actually a functional of the associated 
branching process, n > 0), since the random variable Z„ represents the number of 
individuals at the n-th generation in the forest. The expression of this law was first 
obtained by Dwass |4j. 
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Theorem 1.1 (Dwass, 69). Let Z = {Zn) be a critical or subcritical branching process. 
Let Pfc be its law when it starts from Zq = k > 1 and denote by v its progeny law. Let O 
be the total progeny of Z, that is O = Yln>o ^n- Then for any n > k, 

(1.1) Fk{0 = n) = -iy*^{n-k), 

n 

where z/*" is the n-th iteration of the convolution product of the probability v by itself. 
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Figure 1. A forest labeled according to the breath first search order and 
the associated Lukasiewicz-Harris coding path. 

More generally, whenever a functional of the branching forest admits a nice expression in 
terms of the Lukasiewicz-Harris coding path, we may expect to obtain an explicit form of 
its law. For instance, the law of the number of individuals with a given degree in the k 
first trees can be obtained in this way. We refer to Proposition 1.6 in [BJ where the law of 
the number of leaves, first obtained in (TU], is derived from the Lukasiewicz-Harris coding. 
The aim of this paper is to extend the coding that is recalled above to any irreducible, 
critical or subcritical multitype branching forest, and to derive the joint law of the total 
progeny of each type. Up to now, most of the results in this direction concern non irre- 
ducible, 2-type branching processes, as we will recall after the following definitions and 
notations. 

We set Z+ = {0, 1,2,...} and N = {1, 2, . . . }, and in all the sequel of this paper, d 
will be an integer such that d>2. On a probability space (f2, Q, P), we define a d-type 
branching process Z := {izi'\...,zL'^),n > 0}, Z^|_ valued Markov chain with 

transition probabilities: 

P(Z„+i = {ki, . . . , fc,) I Z„ = in, r^)) = * ■ ■ ■ * u^/^ih, ...,kd), 

where Ui are distributions on and z/"' is the r-th iteration of the convolution product of 
Ui by itself, with z/*° = 6o. For r = (ri, . . . , r^) G Z^^, we will denote by Pj. the probability 
law P( ■ I Zq = r). The vector u = {ui , . . . , Uff) will be called the progeny distribution of 
Z. According to this process, each individual of type i gives birth to a random number of 
children with law Ui, independently of the other individuals of its generation. The integer 
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valued random variable Zn is the total number of individuals of type i, at generation n. 
For i,j e {1, . . . ,d}, let us define the rate 

that corresponds to the mean number of children of type j, given by an individual of type 
I and let 

M := (mij)i<ij<rf 

be the mean matrix of Z. Suppose that the extinction time T is a.s. finite, that is 

(1.2) T = inf {n : Z„ = 0} < cx) , a.s. 

Then let Oj, be the total number of offspring of type i, i.e. the total number of individuals 
of type i which are born up to time T: 

n=0 n>0 

The vector (Oi, . . . , O^) will be called the total progeny of the multitype branching process. 
In the case where d = 2 and when mi2 > and < mn < 1 but m22 = ^21 = 0, it may 
be derived from Theorem 1 (n) in that the distribution of the total progeny of Z is 
given by 

(1.3) P( o)(0i = ni,02 = ^2) = — i^r'('^i - ^^1,^2) , l<ri<ni. 

Hi 

When mi2 > and < mii,m22 < 1 but 77221 = 0, after some elementary computation, 
combining the identities in (11. ip and (11.31) . we obtain that for n2 > 1, 

(1.4) P(.,0)(Ol = 721,02 = 772) = ^ Vj-Z^rn^l " ^ , J>2"^ (0, ^2 - j) • 

721772 ^ 

Note that (II. 3p and (II. 4p concern only the reducible case, when d = 2 and T < 00, a.s. 
As far as we know, those are the only situations where the law of the total progeny of 
multitype branching processes is known explicitly. 

Recall that if M is irreducible, then according to Perron- Frobenius Theorem, it admits 
a unique eigenvalue p which is simple, positive and with maximal modulus. In this case, 
we will also say that Z is irreducible. If moreover, Z is non-degenerated, that is, if 
individuals have exactly one offspring with probability different from 1, then extinction, 
that is (II. 2p . holds if and only if p < 1, see [6], [H] and Chapter V of [Ij. If p = 1, we 
say that Z is critical and if p < 1, we say that Z is subcritical. The results of this paper 
will be concerned by the case where Z is irreducible, non-degenerated, and critical or 
subcritical so that (II. 2p holds, that is the multitype branching process Z becomes extinct 
with probability 1. 

The next result gives the joint law of the total progeny together with the number of 
connected components of each type in the forest, which can also be defined as the total 
number of individuals of type j, whose parent is of type i, i ^ j, up to time T. Let us 
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denote by this random variable. Then Theorem 11.11 and identity fll.4p are extended 
as follows: 

Theorem 1.2. Assume that the d-type branching process Z is irreducible, non- degenerated 
and critical or subcritical. For i = 1, . . . , d, let Oi be the total offspring of type i, up to 
the extinction time T and for i ^ j , let Aij be the total number of individuals of type j , 
whose parent is of type i, up to time T. 

Then for all integers ri, rii, kij, i,j = 1, . . . ,d, such that ri > 0, ri + ■ • • + > 1, 
kij >0,foriy^j, -kjj = rj + Yji^j ^ij, and Ui > -ku, 

Pr(Oi = m, . . . ,0d = Ud, Aij = kij, I < i,j < d,i ^ j) 

= ~ I I '{kii, • • • , fcj j-i, rii + kii, ki i^i, . . . , kid) ■, 

1=1 

where r = (ri, . . . ,rd), z/*° = 6q, fii = UiM 1 and K is the matrix {—kij)i<ij<d to which 
we removed the line i and the column i, for all i such that Ui = 0. 

As in the single type case, it is theoretically possible to obtain the law of the total 
progeny (Oi, . . . , Od) of a critical or subcritical multitype branching process by applying 
the multivariate Lagrange inversion formula, to its generating function. An attempt in 
this direction has been done in [5], but the results are not directly explicit in terms of 
the progeny distribution. Moreover, we emphasize that the random variables Aij are not 
functionals of the multitype branching process Z. So, their formal definition and the 
computation of their law require a more complete information provided by the forest. 

Our proof of Theorem 11.21 uses a bijection, displayed in Theorem 12. 7^ between multitype 
forests and a particular set of multidimensional, integer valued sequences. A consequence 
of this result is that any critical or subcritical irreducible multitype branching forest is 
encoded by d independent, (i-dimensional random walks, see Theorem 13.11 Then, in a 
similar way to the single type case, the total progeny, jointly with the number of subtrees 
of each type in the forest, is expressed as the first passage time of this multivariate process 
in some domain. An extension of the Ballot Theorem that can be considered as one of 
the main result of this paper, see Theorem 13. 4[ allows us to conclude, as in the single 
type case. We emphasize that this coding should allow us to compute the law of other 
functionals of branching forests such as the number of vertices of a given degree. 

This paper is organized as follows. Section |5] is devoted to deterministic multitype 
forests. In Subsection 12. H we define the space of these forests. Then Subsection 12.21 is 
devoted to the definition of the space of the coding sequences and to obtain the bijection 
between this space and the space of multitype forests. This result is stated in Theorem 
12.71 Then in Section |3l we define the probability space of multitype branching forests, 
we display their multitype Lukasiewicz-Harris coding in Theorem 13.11 and we prove its 
application to the total progeny that is stated in Theorem 11.21 This result requires a 
multivariate extension of the Ballot Theorem, see Theorem 13. 4[ whose proof bears on the 
crucial combinatorial Lemma 13.31 The latter is proved in Section |H 
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2. MULTITYPE FORESTS 

2.1. The space of multitype forests. A forest, is a directed planar graph f C v x v, 
with a finite or infinite set of vertices v = v(f), such that the outer degree of each vertex 
is equal to or 1 and whose connected components, which are called the trees, are finite. 
A forest consisting of a single connected component is also called a tree. In a tree t, the 
only vertex with outer degree equal to is called the root of t. It will be denoted by r(t). 
The roots of the connected components of a forest f are called the roots of f . For two 
vertices u and f of a forest f , if {u, v) is a directed edge of f , then we say that u is a child 
of V, or that v is the parent of u. 

A label will be given to each vertex of a forest f in the following way. We first give 
an order to the trees of f and denote them by ti(f), t2(f), . . . , tfc(f), . . . (we will usually 
write ti,t2, . . . ,tfc, . . . if no confusion is possible). Then each tree is labelled according 
to the breadth first search algorithm: we read the tree from its root to its last generation 
by running along each generation from the left to the right. This definition should be 
obvious from the example of Figure [TJ If a forest f contains at least i vertices, then the 
2-th vertex of f is denoted by Wj(f). When no confusion is possible, we will simply denote 
the i-th vertex by -Uj. The set of plane forests will be denoted by J-". The elements of J-" 
will simply be called forests. 

Recall that d is an integer such that d > 2. To each forest f G J-", we associate an 
application Cf : v(f) — t- {1, 2, . . . ,d} such that if Ui, tti+i, . . . , Ui+j G v(f) have the same 
parent, then Cf(uj) < Cf(tti+i) < ■ ■ ■ < Cf(Mj+j). For v G v(f), the integer Cf{v) is called 
the type (or the color) of v. The couple (f , Cf ) is called a d-type forest. When no confusion 
is possible, we will simple write f . The set of d-type forests will be denoted by J-^. A 
2-type forest is represented on Figure [2] below. 




Figure 2. A two type forest labeled according to the breath first search 
order. Vertices of type 1 (resp. 2) are represented in white (resp. black). 
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A subtree of type 2 G {1, . . . , d} of a rf-type forest (f , Cf) G J-^ is a maximal connected 
subgraph of (f , Cf ) whose all vertices are of type i. Formally, t is a subtree of type i of 
(f, Cf), if it is a connected subgraph whose all vertices are of type i and such that either 
r(t) has no parent or the type of its parent is different from i. Moreover, if the parent of 
a vertex v G v(t)^ belongs to v(t), then Cf{v) ^ i. Subtrees of type i of (f , Cf) are ranked 
according to the order of their roots in f and are denoted by , t^*\ . . . , t[^'\ . . . . The 
forest f^*) := {t^^^tg^ . . . jt^*'', . . . } is called the subforest of type i of (f,Cf). It may be 
considered as an element of J-". We denote by Ui \ \ . ■ ■ the elements of v(f'^*)), ranked 
in the breadth first search order of f*^*). The subforests of type 1 and 2 of a 2-type forest 
are represented in Figure El 
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Figure 3. Labeled subforests f^^^ and f^^^ associated with Figure |5J Beside 
each vertex, the first number corresponds to its rank in f*^*\ i = 1,2 and 
the second one is its rank in the original forest f . 



To any forest (f , Cf ) G J-^, we associate the reduced forest, denoted by (fr,Cf^) G J-'d, 
which is the forest of J-^ obtained by aggregating all the vertices of each subtree of (f , Cf ) 
with a given type, in a single vertex with the same type, and preserving an edge between 
each pair of connected subtrees. An example is given in Figure |H 

2.2. Coding multitype forests. For a forest (f , Cf ) G J-^ and u G v(f), when no con- 
fusion is possible, we denote by Pi{u) the number of children of type i of u. For each 
i G {1, . . . , fi}, let > be the number of vertices in the subforest f*^*) of (f , Cf). Then let 
us define the rf-dimensional chain ' ), with length and whose values 

belong to the set Z'^, by Xq^ = and if > 1, 

(2.5) x%,-xli^ =pj{v!^^, ifz^jand x^^^^ = - 1 , 0<n<n,-l, 
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Figure 4. Reduced forest associated with the example of Figure |21 Beside 
each vertex is the rank of the root of the corresponding subtree in the 
original forest f . 

where we recall that (un'')n>i is the enumeration of the subforest f^*) in its own breadth 
first search order. Note that the chains (x'^f), for i ^ j are nondecreasing whereas (x^*) 
is a downward skip free chain, i.e. x^^j^ — x^* > —1, for < n < rii — 1. The chain (x^*) 
corresponds to the Lukasiewicz-Harris coding walk of the subforest f*^*\ as defined in the 
introduction, see also Section 6.2 in [13j for a proper definition. In particular, if Ui is 
finite, then rii = min{n : x^* = mino<fc<ni ^k^}- These properties of the chains x^^^ lead us 
to the following definition. 

Definition 2.1. Let Sd be the set of -valued sequences, x = {x^^\ x^^\ . . . , x^'^^), such 
that for all i = 1,. . . ,d, z*-*-* = (x*'^, . . . is a TJ^-valued sequence defined on some 

interval of integers, {0, . . . , Ui}, < rij < oo, which satisfies Xg'' = and if Ui > 1 then 

{i) fori 7^ j, the sequence (a;^-')o<„<ni is nondecreasing, 
(ii) for all i, x^ll^i — x^* > —1, < n < Ui — 1. 

We will also denote sequences in by [x^/l'' , < k < Ui, i,j = 1,. . . ,d). The vector 
n = (ni, . . . ,77,^) G if^, where Z+ = Z+ U {+00} will be called the length of x. 

Relation (12. 5p defines an application from the set J-^ to the set Sd- Let us denote by 
this application, that is 

(2.6) ^-.J^d ^ Sd 

(f,Cf) ^ ^iii,Cf))=X. 

For X G Sd, set ki = — info<n<ni and define the first passage time process of the chain 
(x^*) as follows: 

(2.7) T^'^ = min{n > : x^f = -k} , < k < h , 

where rj:^ = 00, if ki = 00. If x is the image by \I' of a forest (f , Cf) G J^d, i-e. x = ^^((f , Cf)), 
then ki is the (finite or infinite) number of trees in the subforest f^*) and for k < 00, the 
time r^*"* is the total number of vertices which are contained in the k first trees of f'-*-', i.e. 
t[*\t2'', . . . , t^!'^ This fact is well known and easily follows from the Lukasiewicz-Harris 
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coding of the single type forest f'^*^ see the introduction and Lemma 6.3 in [13]. Then for 
i,j = 1, . . . ,d, define the integer valued sequence 

(2.8) 4^ =a;''^ (rf ), < A: < A;^ . 

If X = \E'((f , Cf)), then we may check that when i ^ j, is the number of subtrees of type 
j whose root is the child of a vertex in t[*\ \ . . . , t^*''. Or equivalently, it is the number 
of vertices of type j whose parent is a vertex of t^*'', t2'*, . . . , t^*"*. For each i = 1, . . . , d, we 

set 

x(*) = . . . , x^'^) and x= (xW,x(2),...,x('=')). 

Clearly for % ^ j, the sequence {x\-')o<k<ki is increasing and = —k, for all z = 1, . . . , d 
and < k < ki, so that x & Sd and recalling the definition of the reduced forest, (fr, Cf^), 
see the end of Section l?!T| we may check that: 

(2.9) $((fr,CfJ) =x. 

For a forest (f , Cf) G J-^ with trees ti, t2, . . . , we will denote by C(f^cf) the sequence of types 
of the roots of ti, t2, . . . , i.e. 

C{f,cf) := (cf(r(ti)),Cf(r(t2)),...). 

Note that C(f cf) G Ui<r<oo{li • • • i f^}*^ and that C(f^(.f) = C(fy,cf^)- When no confusion is 
possible, C(f^cf) will simply be denoted by c = (ci, C2, . . . ) and we will call it the root type 
sequence of the forest. 

Then before we state the general result on the coding of multitype forests in Theorem 
12. 7[ we first need to show that the sequences {x^'^)i-Lj together with c = (ci, C2, . . . ) allow 
us to encode the reduced forest (fr,Cf^), i.e. this forest can be reconstructed from (x,c). 
This claim is stated in Lemma [2.51 below. In order to prove it, we first need to describe 
the set of sequences which encode reduced forests and to state the preliminary Lemma 
12.21 regarding these sequences. 

Recall that Z+ = Z+ U {+C)o} and let us define the following (non total) order in Z^: 

for two elements q = (gi, . . . , g^) and q' = {q[, . . . , q'^ of Z^ we write q < q' if gj < q[ 
for all 2 = 1 , . . . , (i. Moreover we write q < q' if q < q' and if there is z = 1 , . . . , (i such 
that qi < q[. For an element x = {x\\ < k < qi, i,j = 1, . . . ,rf) of with length 
q = {qi, . . . , qa), and r = (ri, . . . , r^) G Z^j., we say that the system of equations (r, x) 
admits a solution if there exists s = {si, . . . , Sd) G Z^, such that s < q and 

d 

(2.10) r, + J]x^'^(s,) = 0, j = l,...,d. 

Lemma 2.2. Let x & Sd and r = (ri, . . . , r^) G Z'|. Assume that the system (r, x) admits 
a solution, then 

(i) there exists a unique solution n = (rii, . . . ,nd) of the system (r,x) such that if 
n' is any solution of (r,x), then n < n'. Moreover we have = min{n : x^* = 
mino<A;<ni x^}: for all i = 1, . . . , d. A solution such as n will be called the smallest 
solution of the system (r,x). 
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(m) Let r' G be such that r' < r. Then the system (r', x) admits a solution. Let 
us denote its smallest solution by n' . Then the system (r — t',x), where x^'^{k) = 
x'^'^{n[ + k) — x^'^{n'j), Q <k <ni — n[, admits a solution, and its smallest solution 
is n — n' . 

A proof of this lemma is given in Section |H For r = (ri, . . . , r^^) G Z*^, with r = ri + ■ ■ ■ + 
r^; > 1, that is r > 0, we define, 

Q = {c G {1, . . . , 4" : Card {1 < j < r : 9 = z} = r^, z = 1, . . . , d} . 

We emphasize that the root type sequence of a forest (f , Cf ) G J^d with r = ri + ■ ■ ■ + 
trees amongst which exactly trees have a root of type i is an element c G C^. Now we 
define the subsets of forests and reduced forests whose root type sequence is in and 
that contain at least one vertex of each type. 

Definition 2.3. Let r = (ri, . . . , r^) G Z'^, such that r > 0. 

{i) We denote by J^^, the subset of of forests (f, Cf) with ri + ■ ■ ■ + rd trees, which 

contain at least one vertex of each type, and such that C(f^cf) ^ C^d- 
(a) We denote by the subset of J-'^, of reduced forests, more specifically, (f , Cf ) G -F^ 
if (f , Cf) G J-'^ and if for each i, vertices of type i = 1, . . . , d in have no child 
of type i. 

Then we define the sets of coding sequences related to J-"^ and 

Definition 2.4. Let r = (ri, . . . , r^) G such that r > 0. 

(i) We denote by S"^ the subset of Sd of sequences x whose length belongs to N'^ and 
corresponds to the smallest solution of the system (r, x) defined in (12.101) . 

(i) We denote by the subset of S]^ consisting in sequences, such that x^ = —k, for 
all k and i. 

Then we first establish a bijection between the sets and S^^ x C^. Recall the definition 
of m in (I23|l . 

Lemma 2.5. Let r = (ri, . . . , r^) G TJ^, be such that r > 0, then the mapping 

(f,Cf) ^ (*((f,Cf)), 

is a bijection. 

Proof. Let (f , Cf ) G and let ki be the total number of subtrees of type i which are 
contained in (f, Cf). (Note that since (f , Cf ) is a reduced forest, its subtrees are actually 
single vertices.) By definition, k^ > 1, for each i. The fact that, C(f^cf) ^ follows 
from Definition 12.31 (ii). Then let us show that x = \l/((f, Cf)) G S"^. Since (f , Cf ) is a 
reduced forest, then x = x. Besides, from (12. 8p . x has length k = {ki, . . . , kd) and for 
j 7^ i, x^'^{ki) is the number of subtrees of type j whose root is the child of a vertex 

of jt^*'', . . . , t^*''|, i.e. of any subtree of type i in (f, Cf). Hence for j = l,...,d, 
x*'-'(/cj) is the total number of subtrees of type j in (f,Cf), whose root is the child 
of a vertex of type i = 1, . . . , d, i ^ j . Then we obtain kj by adding to Yli^tj ^^'''{^di the 



10 



LOIC CHAUMONT AND RONGLI LIU 



Tj subtrees of type j whose root is one of the roots of ti, . . . , t^, where r = ri + ■ ■ ■ + r^. 
That is kj = rj + Xli^j ^*'"' (^«)- Since moreover from fl2.8p . x^'* = —k, for all < A; < fcj, 
we have proved that k is a solution of the system (r, x). It remains to prove that it is the 
smallest solution. 

Let us first assume that r = (1, 0, . . . , 0), so that (f , Cf) consists in a single tree ti whose 
root has color Cf(r(ti)) = 1. Then we can reconstruct, this tree from the d sequences 
{x^j^\ < k < ki), i = l,...,d hj inverting the procedure defined in (12.51) and this 
reconstruction procedure gives a unique tree. Indeed, by definition of the application 
each sequence {x^^\ < k < ki), is associated to a unique 'marked subforest', say 
f^*\ of type i whose vertices kept the memory of their progeny. More specifically, for 
k = 1, . . . , ki, the increment x^*"* — z^*!^ gives the progeny of the A;-th vertex of the subforest 

f*^*\ This connection between marked subforests f^*^ and sequences {x^^\ < k < ki) is 
illustrated on Figure [51 Now let k' < k be the smallest solution of the system (r, x). Let 
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X^^>: AxY' = (-1, 2, 1), A4'' = (-2, 2, 1), Ax]l' = (-3, 2, 1) 
A.x!," = (-4,2,1)- 
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c(2): Axf = (1,-1,0), Axf' = (1,-2,0), Ax.f' = (1,-3,0), 
Axf^ = (1, -4, 2), Axf = (1, -5, 2). 



f(3)^ 



% / 
\ / 



(;(3): Aif = (0, 1, -1), Axf = (1, 1, -2), Axf = (2, 2, -3). 



Figure 5. A three types (reduced) forest (f, Cf), its three marked sub- 
forests and its coding sequences x^^\ x^'^\ x^^\ Here we have 



q = {qij - - - ^Qd) < k' and suppose that we have been able to perform the reconstruction 
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procedure until q, that is from the sequences {x^ ,0 ^ k < qi), i = 1, . . . ,d. Then since q 
is not a solution of (r, x), we see from what has been proved just above that the tree that 
is obtained is 'not complete'. That is, at least one of its leaves (say of type j) is marked, 
so that this leaf should still get children whose types and numbers are given by the next 
jump Xg^^^i — x^g^\ for qj < kj, according to the reconstruction procedure. Thus, doing so, 

we necessarily end up with a tree from the sequences (x^*'', < k < k'^) , i = 1 , . . . , d, and 
this tree is complete, that is none of its leaves is marked. Then since the reconstruction 
procedure obtained by inverting (12.51) . gives a unique tree, we necessarily have k' = k. 

Then let r = (ri, . . . , r^) G Z"^. Assume with no loss of generality that the root of 
the first tree ti of (f , Cf ) has color 1. Let kj be the number of subtrees of type i in ti. 
From Lemma [2. 2^ the system (r^, x), where r^ := (1,0,..., 0), admits a smallest solution. 
Moreover from the reconstruction procedure which is described above, this solution is 
k^ = {kl,...,k]). Suppose now with no loss of generality that the second tree, t2 in 
(f , Cf ) has color 2. Let kf be the number of subtrees of type i in t2. Then from the 
same arguments as for the reconstruction of the first tree, t2 may be reconstructed from 
the system (r^, y), where r^ := (0, 1, 0, . . . , 0) and y^'^{k) = x^'^{kl + k) — x^'^{k}), k > 0. 
Moreover {I'^^y) admits k^ = (A;^, . . . , /c^) as a smallest solution. Then from part {ii) of 
Lemma [221 k^ +k^ is the smallest solution of the system (r^ +r^, x). So we have proved the 
result for the forest consisting in the trees ti and t2. Then by iterating these arguments 
for each tree of (f , Cf), we obtain that x G S"^. 

Conversely, let c = (ci, . . . , Cr) G C^, x G and let k = (A;i, . . . , kd) be the small- 
est solution of the system (r,x). Then let us show that there is a forest (f , Cf ) G 
such that \E'((f, Cf)) = x and C(f^cf) = c. Assume, without loss of generality that ci = 1. 
From Lemma 12.21 {ii), there is a smallest solution, say k^ = (k\, . . . , fc^), to the system 
(r^, x), where r^ := (1, 0, . . . , 0). Then we may reconstruct a unique forest (ti, Ct J G ^\ 
(consisting in a single tree) such that \E'((ti,CtJ) = ^ k < k], i,j = l,...,d) 

and C(tj,ctJ = 1 by inverting the procedure that is described in 02. 5p . Assume for in- 
stance that C2 = 2 and set r^ := (1, 1, 0, . . . , 0) < r, then from Lemma 12.21 {ii), there is 
a smallest solution, say k^, to the system (r^,x). Moreover, k^ = k^ — k^ is the smallest 
solution of the system (r^ — T^,y), where y^'^{k) = x^'^{kl + k) — x'^'^{k}), k > 0. Then 
as before, we can reconstruct a unique tree (t2,Ct2) such that \E'((t2, Ctj)) = y and such 
that the forest f = {ti,t2} satisfies \E'((f,Cf)) = (a;^'"',0 < k < kf, i,j = l,...,d) and 
^(fc-) ~ (1) 1) 0) • • • ) 0). Then iterating these arguments, we may reconstruct a unique 
forest (f , Cf) G J-"^ such that \l/((f , Cf)) = x and C(f^cf) = c. □ 

Let X & Sd with length n = {ui, . . . ,nd) and recall from f l2.8p . the definition of the 
associated sequence x, with length k = {ki, . . . , k^), such that ki = — mino<„<ni x^^. 

Lemma 2.6. Let r G TJ^, such that r > and x G Sd, with length n G N'* and set 
ki = — mino<n<ni x^i^ , i = 1, . . . ,d. If n is the smallest solution of the system (r, x) {i.e. 
X G S"^), then k = {ki, . . . , kd) is the smallest solution of the system (r, x). Conversely, if 
rii = tI:\ for all i = 1, . . . , d and ifk is the smallest solution of (r, x), {i.e. x G S'^), then 
n is the smallest solution of (r, x). 
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Proof. Assume that n is the smallest solution of the system (r, x). Then from part (i) 
of Lemma fI72\ rii = r^*\ hence k is a solution of (r, x). Let k' < k be such that 

d 

1=1 

Then by definition of x there is n' < n such that n[ = t^) and 

d 

rj + ^x''^{n',)=Q, j = l,...,d. 

i=l 

So n' = n and hence k' = k. 

The converse is proved in the same way. Suppose that rij = t^*\ i = 1, . . . , d, and that 
k is the smallest solution of (r,a;). Then clearly, n is a solution of (r,a;). Let n' be the 
smallest solution of (r,x). Then from Lemma [2. 2[ there is k' such that n- = Tu\ hence k' 

i 

is a solution of (r, x). This implies that k < k', so that n < n', hence n = n'. □ 

Now we extend the application $ defined in Lemma [2.51 to the set J-"^. Here is the main 
result of this section, that can be considered as an extension of Proposition 1.1 in [3]. 

Theorem 2.7. Let r = {ri, . . . ,rj) G Z^, be such that r > 0, then the mapping 

(f,Cf) ^ (vl/((f,Cf)),C(f,,,)) 

is a bijection. 

Proof. Let us first check that for any (f , Cf ) G J-"^, we have $((f, Cf)) G 5"^ x C^. By 
definition 12. 3[ (i), C(f cj.) G C^. Now set x = \l/((f, Cf)) and let (fr,Cf^) G be the forest, 
(f , Cf ) once reduced. Then from (12. 9p and Lemma 12. 5[ this reduced forest is encoded by 
{x, Cf). Let k = (A;i, . . . , k^) be the number of subtrees of type i in this forest (this is 
actually the number of vertices of type i), then k is the length of x and it is the smallest 
solution of (r,x), i.e. x G S"^. Moreover n = {ni, . . . ,nd), where = rj: ^ is the length 
of X, and from Lemma [2.6[ it is the smallest solution of (r,x). So, we have proved that 
^((f,Cf))G55. 

Conversely let {x, c) G x C^. From Lemma [2.51 to {x, c), we may associate a unique 
forest (fr,Cf^) G Then let ki be the number of vertices of type i in this forest. For 
1 < A; < A;j, let be the A;-th vertex of type i in the breadth first search order of (fr, Cf^). 
Then in (fr, Cf^), we replace the vertex u^^^ by the subtree of type i which is encoded by the 
Lukasiewicz-Harris path (x*'*(r^*2]^ + /) + A; — 1, 0</< r^*^ — t^*2i). We know about the 
progeny of each vertex of this subtree, thanks to the chains {x'^'-' (t^Ii + I) — x^'^ {Tj:^-^) , < 

/ < r^*'' —rj^^li), so that we can graft at the proper place, on this subtree, all the correspond- 
ing subtrees of the other types which have been constructed from the same procedure. 
Proceeding this way, we construct a unique forest (f , Cf ) G J-"^ and we easily check that 

^((f,Cf))=X. □ 
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3. MULTITYPE BRANCHING TREES AND FORESTS 

Let Ui, i = 1, . . . , d he distributions on Z^, such that u = [ui, . . . ,i/d) is the progeny 
law of an irreducible, critical or subcritical, non-degenerated branching process, as de- 
fined in Section [U Assume that we can define on the reference probability space {Q, Q, P) 
introduced in Section [1], a family (Pc)cG{i,...,rf}°° of probability measures and an infinite 
sequence F = (Tfc)/c>i of independent random trees, such that for each c = (ci, C2, . . . ) G 
{1, . . . , d}'^ and A; > 1, under P^, is a branching tree, with progeny law u, whose root 
has type r(Tfc) = c^. In particular, for any random time a : {Q,Q) — ?■ N U {+00}, the 
sequence {Ti, . . . , T^} is an element of J-'^. The infinite sequence F will be called a d-type 
branching forest with progeny law ly. 

Let us denote by F*^*-' the subforest of type i of F, as it is defined in subsection 12. II From 
the properties of u, it follows that for each i = 1, . . . ,d, the subforest F*^*^ is a.s. infinite, 
so that we may define a Z'^ valued infinite random sequence 
i = 1, . . . , d, in the same way as in (12. 5p . that is X^^ = and 

(3.11) JCii-JQ^=p,(u^|i), ifz^jand - = p,(uj|i) - 1 , n>0, 
where (un'')„>i is the enumeration of F*-*) in its breadth first search order. 

Theorem 3.1. Let F be a d-type branching forest with progeny law u. 

1. Then for any c = (ci, C2 . . . ) G {1, . . . , d}°°, under P^, the chains 

(3.12) X^'^ = {X''\ . . . , X'''^) , i = l,...,d 
are independent random walks with step distribution 

Pc{XP = (gi, . . . , qa)) = Vi{qi, gi_i, Qi + 1, g^+i, ...,qd). 

In particular, their laws do not depend on c. For each i = 1,. . . ,d, X*'* is a 
downward skip free random walk such that lim inf „^_|.oo AT^'* = —00, a.s. and for 
j 7^ i, X*'-^ is a renewal process. 

2. For all integer r > 1, almost surely there is t E Z'j^, with r = ri + ■ ■ ■ + r^ and 
such that there is a smallest solution n to the system {t,X). 

3. Conversely, let Y be a copy of X and c = (ci, C2 . . . ) G {1, . . . , d}°° . Then to Y 
and c, we may associate a unique d-type branching forest, with progeny law 1/ and 
root type sequence c, whose coding random walk is Y . 

Proof. Part 1. just follows from the construction (13. lip of X. Since the order on 
the subforests F^*^ does not depend on the particular topology of F, from the branching 
property, it is clear that the chains i = l,...,d are independent random walks. 

Then the expression of the law of X*^*) is a direct consequence of (13.111) . Recall that X*'* is 
the Lukasiewicz- Harris path of the subforest F^^\ see section 6.2 of [13j. Moreover, since 
from the properties of u, each subforest F^*) is a.s. infinite, the random walk X*'* satis- 
fies lim inf„_!.oo X*'* = —00, a.s. The fact that X*'-' , for z 7^ j is a renewal process is obvious. 

Then part 2 is a direct consequence of the construction of X and Theorem 12.71 Let 
r > 1 and first assume that the finite forest {Ti, . . . , T,,} consisting in the r first trees of 
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F contains at least one vertex of each type. Let r be the unique element of Z'^ such that 
{Ti, . . . ,Tr} G J-"^. Then, by coding the forest {Ti, . . . , T^.} and by applying Theorem 
\2.7\ we obtain that there is n G N"^ which is the smallest solution of the system {r,X). 
Now if for instance n = 1, . . . , d — 1 types are missing in the r first trees of F, then we 
can apply the same arguments by replacing dhy d — n in Theorem 12.71 

Then part 3. is a consequence of part 2. For each r > 1, we may associate a unique 
forest to Y with r trees. Since r can be arbitrarily large, the result is proved. □ 

In the same spirit as in [9j, the Lukasiewicz-Harris type coding that is displayed in 
Theorem 13.11 might be used to obtain invariance principles, for any functional that can 
be encoded simply enough. Besides this result should provide a way to obtain a proper 
definition of continuous multitype branching trees and forests. Actually, it is natural to 
think that the latter objects are coded by d independent, d- dimensional Levy processes, 
with d — 1 increasing coordinates and a spectrally positive coordinate. 

Now we are going to apply our coding of multitype branching forests to the law of their 
total progeny and give a proof of Theorem 11.21 To that aim, we first need to establish 
the crucial combinatorial Lemma 13.31 Let E be Z+ or a finite integer interval of the type 
{0, 1, . . . ,m}, with m > 1 and let g : E ^ Z'^, be any application such that g{0) = 0. 
For n & E such that n > 1, the n-cyclical permutations of g are the n applications gq^n, 
g = 0, . . . , 77, — 1 which are defined on E by: 



(3.13) g^h) ^ = ^ 



g{q + h) — g{q) if < h < n — q , 

g{h — {n — q)) + g{n) — g{q) if n — q < h < n 
g{h) if h> n . 



Note that go^n = 9- The transformation g g^ ,^ consists in inverting the parts {g{h), < 
h < q} and {g{h), q < h < n} in such a way that the new application, gg^n, has the same 
values as at and n, i.e. 5'g,n(0) = and g^nin) = g{n). 

Let X G Sd, with finite length n = (ni, . . . , n^^) G N'^ and recall the notation x^^^ = 
(x*'^, . . . , x*''^) from Definition 12.11 Then we define the n-cyclical permutations of x by 

(3.14) Xq^n := (4ilni> • • • > ^i^nj ' ^ = (^1' ■ ■ ■ , Qd) < ^ - U, 

where we have set Id = (1, 1, . . . , 1). Each sequence Xq_n will simply be called a cyclical 
permutation of x. Note that there are nin2 . . . Ud, cyclical permutations of x. Let r = 
(ri, . . . , r^) G 'Z'il be such that r > and assume that n is a solution of the system (r, x). 
Then note that n is also a solution of the system (r, x^^^), for all q < n — 1^^, that is. 



i=l 



This remark raises the question of the number of cyclical permutations n of x, such 
that n is the smallest solution of the system (r,Xq^n)- 
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Definition 3.2. Letx G Sd, with finite lengthn = (rii, . . . ,nd) G N'^. Leti = (ri, . . . ,rd) & 
be such that r > and assume that n is a solution of the system (r, x) . For q < n — 1^, 
we say that good n-cyclical permutation of x with respect to r, ifn is the smallest 

solution of the system (r, Xq^n); that is Xq^n ^ S^^. When no confusion is possible, we will 
simply say that good cyclical permutation ofx. 

Here is our crucial combinatorial cyclic lemma. 

Lemma 3.3. Let x G Sd, with length n = (ni, . . . ,77,^) G N'' and let r = (ri, . . . , r^) G 
be such that r > 0. Assume that n is a solution of the system (r, x) such that 
x^'^{ni) 7^ 0, for all i = 1, . . . ,d. Then the number of good cyclical permutations of x 
is det((-x''^(ni))i<ij<d). 

This lemma will be proved in Section HI It is used for the proof of the following extension 
of the Ballot Theorem. 

Tiieorem 3.4 (Multivariate Ballot Theorem). Let Y = (Y^^\ . . . ,Y^'^^) be a stochastic 
process defined on {n,g,P), with = (y*'J(/i),l < j < d, h > 0) , i = 1, . . . , d and 
Yq = 0. We assume that the coordinates Y^'^ , for i ^ j are Z+ valued, nondecreasing and 
that the coordinates F*'* are Z valued and downward skip free. Fix n = (ni, . . . , n^) G N'^, 
then we assume further that the process Y is n- cyclically exchangeable, that is for any 
q = (gi, . . . , gd) G Z'l^ such that q < n — 1^, 

V ) V 

q,n — 5 



where Fq^n is defined as in f l3.14p for deterministic functions . Then for anyi = (ri, . . . , r^) G 
T<i oni^h +hn+ T ^ Q and kij, 1 < i,j < d, such that kij G Z+, for i ^ j and —kjj = 



P{Y^f = kij, 1 < i, j < d and n is the smallest solution of (r, Y)) 

(3.15) = _Mzh2lp(^Y^^ = k,,, = hj, ^^J,l<^,J<d). 

Proof. If P{Y^'^ = kij,i ^ j) = 0, then the result is clearly true. Suppose that it 
is not the case and let y = {y^^\ . . . ,2/*-'^'') be a deterministic function such that for all 
^ < hj < d, y^'^ is defined on Z+, (0) = 0, y'^'^iui) = kij and 

(3.16) P{{Y^'\h),0 <h< m)) = {y^'\h),0 <h< m)) > . 

For h = {hi, ...,hd)e Z^, we set Y{h) := (rW(/ii),. . . ,yW(/id)) and for < q < n-l^, 

using the notation of (I3.14p . we set 2/q,n(h) := {ycn]ni{hi), . . . , yqf,na{hd)) ■ Let us consider 
the set 

Ey,n = {(?/q,n(h), < h < u) : < q < u - 1 } , 

of n-cyclical permutations of y over the multidimensional interval [0, n]. Then Caid{Ey^^) = 
nin2 ■ ■ - Ud and since Y = {Y^^\ . . . , Y^*^^) is a cyclically exchangeable chain, the law of 
(y(h), < h < n), conditionally to the set {(y(h), < h < n) G Ey^^} is the uniform 
law in the set Ey^^. Moreover, assume that ku ^ for all i = 1, . . . ,d, then conditionally 
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to the set {(y(h), < h < n) G Ey^^}, from Lemma 13.3^ the number of good cychcal 
permutations of {Y{h.), < h < n) is det{—kij). Therefore, 

P(y^'"' = kij, I < i, j < d and n is the 

smallest solution of (r, Y) I (y(h), < h < n) G ^„„) = ^ ^ . 

nin2 ...Ud 

Then we obtain the result by summing the identity 
P{Y^f = kij, l<i,j <d and 

n is the smallest solution of (r, Y), {Y{h), < h < n) G Ey^^) 



det{- kij] 



-F((y(h),0<h<n)G^,,n), 



nin2 ...Ud 

over all functions y satisfying f l3.16p and with different sets Ey^^ of cyclical permutations. 
Finally, if ka = 0, for some i = 1, . . . , d, then since > 1, we can see that both members 
of identity fl3.15p are equal to 0. □ 

Proof of Theorem 11.21 Let r, n and kij be as in the statement. Let F be a d-type 
branching forest with progeny law u, as defined at the beginning of this section and such 
that the r first trees have root type sequence (ci,...,Cr) G C^. Let X be the coding 
random walk of F, as defined in (13. lip . Recall the notation of Theorem II. 2[ then from 
the coding of Subsection 12.21 and Theorem 13. 1[ we may check that 

(3.17) Pr(Oi = m,. . . ,0d = nd,Aij = kij, l<i,j <d,i^ j) 

= Pc(-^^f = kij, 1 < i, j < d and n is the smallest solution of (r,X)) . 

Assume first that n G N'^, then since X is clearly cyclically exchangeable in the sense of 
Theorem 13. 4[ we obtain by applying this theorem, 

Pc(X^f = kij, 1 < i, j < d and n is the smallest solution of (r, X)) 

= ^"^(-^-^^ P(x^; = ku, x^ = h„i^ j,i<i,j<d). 

niT2,2 ■ ■ - rid 

On the other hand, since from Theorem 13. 1[ the random walks X^^\ i = 1,. . . ,d are 
independent, we have 

Pc(^^f = kij, 1 < i, j < d and n is the smallest solution of (r, X)) 

= ^^^^^tlp{xi;: = k.,xi;^ = h„z^j,i<j<d). 

nin2...nd^^ 

1=1 



Then from the expression of the law of X given in Theorem 13. 1[ we obtain 

VciXli^, = kij, 1 < i, j < d and n is the smallest solution of (r, X)) 

(3.18) = - I I h'j ^ (kii, . . . , ki i^i, rii -\- kii, ki i-^-i, . . . , ki d) , 

and the result is proved in this case. 
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Now with no loss of generality, let us assume that for some < rf' < rf, we have 
ni,...,nd' G N and nd>+i = ■ ■ • = = 0. We point out that from the assumption 
rij > —kjj = Tj + Yli^j kij, in this case we necessarily have rj = and kij = 0, for 
all i = 1, . . . ,d and j = d' + 1, . . . ,d. Then provided we also have kij = 0, for all 
i = d' + 1, . . . ,d and j = 1, . . . ,d, i ^ j, 

p^(X;f = kij, d' + l<i<d, l<j <d) 

d 

— ^{kiii . . . , ki^i—i,Tii + kii, ki^i^i, . . . , ki^d) = 1 • 

i=d'+l 

Define the chain X restricted to Z'^', by X' = {X'^^\ where = {X''\ . . . , X''"^'] 

i = 1, . . . ,d'. Set also n' = (ni, . . . ,nd>) and r' = (ri, . . . , r^/). Then under our assumption 
on the integers kij, the following identity is satisfied, 

{Xl^^ = kij, 1 < i, j < d and n is the smallest solution of (r,X)} 
= {X'^f = kij, 1 < i,j < d' and n' is the smallest solution of (r', X') and 

= o,i<i<d',d' + i<j<d}, 

so that identity fl3.17p can be rewritten as, 

F,{Oi = m,. . . ,0d = Ud, Aij = kij, I <i,j < d,i ^ j) 
= Pc{X'^',^ = kij, 1 < i, j < d' and n' is the smallest solution of (r', X') and 
= 0,l<i<d',d' + l<j<d). 
Moreover, conditionally on the set 

{x:if = h,,, xiii = Q,i<i<d',i<j'<d',d' + i<j<d}, 

the chain X' is cyclically exchangeable, so that we can conclude in the same way as above 
that 

Pi.(Oi = ni,. . . ,0d = rid, Aij = kij, I <i,j < d,i ^ j) 
_ det{{-kij)i<ij<d') 



nin2 ...Ud' 



y ^i ^{kii, . . . , kii_i, rii + kii, ki^i^i, . . . ,ki 



_ det{{-kij)i<ij<d') yj .*n,(j , h h -\ 

nin2...nd' 

1=1 

Finally, if kij ^ 0, for some i = d' + 1, . . . ,d and j = 1, . . . ,d, then the first and the third 
members of the above equality are equal to 0. So the proof is complete. □ 

4. Proof of Lemmas [22] and E3] 

Proof of Lemma 12. 2[ Assume that there is a solution s = {si, . . . , Sd) to the system 
(r, x), that is: rj + X]i=i ^*'"'("^«) ~ 0' J ~ ^,---,d. Let us write this equation in the 
following form: 

+ ^'''(^^) + ^'''(^^■) = ' J = i> ■ ■ ■ > • 
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Then recall that for fixed j, when the ki^s increase, the term Xli^j increases and 

when kj increases, the term x^'^{kj) may decrease only by jumps of amplitude —1. 

For k = {ki, . . . ,kd) = 0, we have rj + ^i^j x''-'^{ki) + x^'^{kj) = rj, j = 1, . . . , d. So 
for the left hand side of the later equation to reach 0, each ki has to be at least r'^^^Vi), 
where r^^^ has been defined in (12. 7p . (In this proof, we found more convenient to use the 
notation T^^\k) for t^*''.) Then either 

r,- + 5^x*'^(r«(r,)) + x^'^'(r(^')(r,)) = , j = l,...,d, 

or all of the terms rj + 'Yliijtj x^'^ (r^^^ri)) + x^'^ {T^^\rj)), j = 1, . . . , d are greater or equal 
than 0, at least one of them being strictly greater than 0. 

Then in the latter case, for rj + Yli^^j x'^'-^{ki) + x^'^{kj) to attain 0, each of the k/s has 
to be at least T^'^\rj + Xli^^j 2;*'-' (T(*)(ri))). This argument can be repeated until all of the 
terms Vj + ^i^j x'^'^{ki) + x^'^{kj) attain 0. More specifically, set v^ ^ = rj and for n > 1, 

^Si = ^. + E^''(^^^^(^"^))' 

and set kj^'' = T^^\vrP). For n > 1, either 

+ 5Z 3^''' (^f^ ) + ^''Kkf ) = , J = 1 , . . . , d , 

or all of the terms Vj + Xli^^j x'^'^{k'f'^) + x^'^{kj^''), j = 1, . . . , d are greater or equal than 
0, at least one of them being strictly greater than 0. In the later case, for all of the 
terms rj + Yli^j x^'^{ki) + x^'^{kj), j = 1, . . . , d to vanish, the index k has to be at least 

= (^k["~^^\ . . . , fc^"^^'')- But since there is a solution s to the equation (r, x), there is 
necessarily a finite index Uq such that k*^"°^ < s and 

rj + Y,x'Hkl"''^) + x"{kf'^) = . 

That is, for all j = k'f"^ = r^^^rj + Ei^j 2;*'J(A;f "^)). Hence k("") is the 

smallest solution of the system (r,x). Moreover by definition, kl"'°'' = min{n : x^* = 
minp^^^^(no) x^/f}. This proves the first part of the lemma. 

Let t' = {r[, . . . , r'^) G be such that r' < r. Then we prove that there is a smallest 
solution to the system (r', x) similarly. More specifically, since there is a smallest solution 
to the equation 

^'j + Yl ^'''(^^) + ^'"'(^j) = r'j-rj<^^ j = l,...,d, 
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then by the same arguments as in the first part, we prove that there is a smallest solution 
to the equation 

r;. + J2 (^0 + ^"(^j) = , 3 = l,...,d. 
Let k and k' be respectively the smallest solutions of (r, x) and (r',x). Then we have, 

d d 

rj - r'. + rr + J2 ^''' (^0 + Yl ^'''^^^ " = ■ 

i=l i=l 

Since r'j + Ylt=i ^''''^ i^'i) — O5 ^^e above equation shows that k — k' is a solution of the 
system (r — r', ot). Moreover, if k" was a strictly smaller solution of (r — r', x), than k — k' 
(i.e. k" < k — k') then we would have 



d d d 

i=l i=l i=l 

= 0, 



so that k' + k" would be a solution of (r, x), strictly smaller than k, which is a contradic- 
tion. □ 



Proof of Lemma 13.31 Recall from Lemma [2.51 that to each forest (f , Cf ) G we can 

associate a system (r,x), with smallest solution k, where x G <S^. We then define the 
Laplacian matrix of each forest (f , Cf ) G by {kij) = {x^'^{ki)). We will also sometimes 
say that (kij) is the Laplacian matrix of the associated coding sequence x. (Actually, in 
the conventional terminology of graph theory, {kij) would rather be the Laplacian matrix 
of a directed graph on the set of vertices {1, . . . , d}, with incidence matrix U = {uij)i<ij<d, 
where Uij = kij, ii i ^ j and = 0, 2 = 1, . . . , (i.) Recall that —kjj = rj + ^i^j kij and 
that for i 7^ j, kij is the total number of vertices in v(f) with type j, whose parent has 
type i and that —ka is the total number of vertices of type i in v(f). 

Definition 4.1. Let i E Z'^ be such that r > 0. 

1. An element (f , Cf) of is said to be a simple forest if for each type i = 1, . . . ,d, at 
most one vertex of type i in v(f) has children {the others are leaves). We assume 
moreover that the types of the roots of simple forests are ranked in the increasing 
order. 

2. For X G <S^, the system {t,x), with smallest solution k = {ki, . . . , kd) is called 
a simple system if for all i = 1, . . . ,d, there is k[ < ki — 1, such that for all 
A; = 0, . . . , fc- and i ^ j, x^'^{k) = and for all k = {k'^ + 1) V ki, . . . , ki and i ^ j , 
x*'^{k) = x''^{ki). 

In other words, simple systems are elements x G iSJ such that for all i = 1, . . . ,d, the 
sequences x^'^ , j = 1, . . . d have at most one positive jump that occurs at the same point 
in {0, ... , ki}. Then we have the following straightforward result. 

Proposition 4.2. A forest of is simple if and only if its associated system {as defined 
in Lemma ^l^ is simple. 
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Definition 4.3. An elementary forest is a forest of J^d that contains exactly one vertex 
of each type. In particular, each elementary forest contains exactly d vertices and is coded 
by the d couples {ji,i), i = 1, . . . ,d, where ji is the type of the parent of the vertex of type 
i. If the vertex of type i is a root, then we set ji = 0. We define the set D of vectors 
{ji, ■ ■ ■ , jd) , ^ ji ^ d such that {ji,i), i = 1, . . . ,d codes an elementary forest. 

Lemma 4.4. Let r = (ri, . . . , r^) G 'Z'^ such that r > and let (r, x) be a simple system 
with smallest solutionk = {ki, . . . , kj) G N'^, then the number of good cyclical permutations 
of x is 

d 

{jl,-,jd)€D i=l 

where for i, j = 1, . . . ,d, kij = x^'^{ki) and where we have set k^i = ri. 

Proof. The proof will be performed by reasoning on forests. Let (r, x) be as in the 
statement. We first make the additional assumption that for all i = 1, . . . ,d, there is j ^ i 
such that kij > 0. Then observe that under this assumption, all the cyclical permutations 
of X are different from each other, since each sequence a;*'-', for i ^ j has exactly one 
(strictly positive) jump. 

Then note that for each good cyclical permutation Xq^k of the system (r, Xq^k) is simple 
itself. Thus from Lemma 1^31 and Proposition 14. 2^ to each good cyclical permutation Xq k of 
X, we may associate a unique simple forest of since, by definition, root type sequences 
of simple forests are ranked in the increasing order. Moreover if q 7^ q' and q, q' < k, 
then from our additional assumption, Xqk 7^ a^q'.k? so that the associated simple forests 
are different. Conversely the simple system (r, y) which is associated to each simple forest 
with Laplacian matrix (kij), through Lemma 12.51 is necessarily obtained from a good 
cyclical permutation of x. Indeed the corresponding sequences y^'^ , for i ^ j, have exactly 
one jump and are such that y'^'^{ki) = kij. So y is nothing but a cyclical permutation of 
X. These arguments prove that the number of good cyclical permutations of x is equal to 
the number of simple forests with Laplacian matrix (kij). 

Then let us prove that the number of simple forests with Laplacian matrix (kij) is 
^(ii jd)&D nf=i First observe that to each simple forest, we can associate a unique 
elementary forest in the following way: the vertex of type j in the elementary forest is 
the parent of the vertex of type i if, in the simple forest the parent of the vertex of type 
i who has children has type j (recall that j = if the vertex of type i is a root). An 
example of an elementary forest associated to a simple forest is given in Figure [6l Then 
let (ji, . . . ,j(i) G D. We easily see that the monomial Y[t=i i^ number of simple 
forests such that for each i, the parent of the vertex of type i who has children has type 
ji. Indeed, there are fcj. j possibilities to choose the vertex of type i who has children. 
In other words, Y[i=i is th^ number of all possible simple forests to which we can 
associate the same elementary forest which is coded hj {ji,i), i = 1, . . . ,d. Then in order 
to obtain the total number of simple forests with Laplacian matrix (kij), it remains to 
perform the summation of these monomials over all the possible elementary forests. So 
we obtained the formula of the statement, under the additional assumption made at the 
beginning of this proof. 
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Figure 6. An elementary forest and its associated simple forest. 



Now assume that for all z = 1, 2, . . . , — 1, there is j ^ i such that kij > and that 
kdj = 0, for all j = 1, . . . , d. Then in any forest with Laplacian matrix (kij), all vertices of 
type d are leaves. Since these leaves are undistinguished, such a simple forest is determined 
only by the relative positions of the vertices of type 1, 2, . . . , c? — 1. So we derive from 
what has just been proved, that the number of simple forests with Laplacian matrix {kij) 
is hd-i := i)e-D' Ilti where D' is the set that is defined in Definition 14.31 and 

where we replaced d hj d — 1. To these hd-^i forests correspond hd-i simple systems to 
which correspond hd-i good cyclical permutations of x. Then in order to obtain all the 
good cyclical permutations of x it remains to consider the kd cyclical permutations of the 
sequence x^'^\ Since the latter are all identical, we have actually kd x hd^i good cyclical 
permutations of x. Then we conclude by noticing that kd x hd-i = jd)&D Y[i=i 

since kd = rd+J2i=^d ^id = Ylt=o,ijtd ^id- The general case where % = 0, for all j = 1, . . . , d, 
for other types i is obtained in the same manner. □ 

Lemma 4.5. For any r G and any integer valued matrix {kij)i<i,j<d, such that kij is 
nonnegative for i ^ j and —kjj = Tj + 'Ylii^j kij, j = ^, ■ ■ ■ ,d, 



where the set D is defined in Definition \4.'d\ and k^^i = ri. 

Proof. The proof is a direct consequence of the matrix tree theorem for directed graphs, 
due to Tutte [13], see Section 3.6, page 470 therein. However, Theorem 3.1 in [12j that 
implies Tutte's theorem is actually easier to apply, since it uses a setting which is closer 
to ours. Let us consider a set {vq, Vi, . . . , Vd} of d + 1 vertices and in the notations of |12] . 
set W = L = {vq}. Then the family J^w,l that is described in Theorem 3.1 of [T^ is in 
bijection with the set of elementary forests, or equivalently with the set D, and identity 
(3.2) in this theorem is exactly det(— A;jj) = J2(n,...,j,)eD 



d 
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Lemma 4.6. Let x E Sd, with length k = {ki, . . . , kd) € N'^- Let r = (ri, . . . , r^) G Z'^ be 
such that r > and assume that k is a solution of the system {t,x). Assume moreover 
that x^'* = —k, for all k = 0, . . . , ki and i = 1, . . . ,d. Then the number of good cyclical 
permutations of x is detl^x"^'^ {ki))i<ij<d) ■ 

Proof. Let x be as in the statement. If for a\\ i = 1, . . . ,d, the sequences a;*'-^, j = 1, . . . d 
have at most one positive jump that occurs at the same point in {0, . . . , ki}, then from 
Definition 14.11 Lemma 14.41 and Lemma 14. 5[ there is a cychcal permutation x' of x such 
that (r, x') is a simple system and the result follows in this case. In general, let us prove 
that there is a simple system, with Laplacian matrix {x'^'^{ki)), and whose number of good 
cyclical permutations is the same as this of x. 

Fix any index m = 1, . . . ,d and assume without loss of generality, that m 7^ 1 , and 

^ ' {km) X ' {km 1) > 0. 

From X, we define a new sequence x E Sd with length k as follows: 

' i:-^'^{k)=x^'^{k), j = 2,...,d, k = l,...,km, 

(4.19) <! x™'i(A;) =x"^'i(A:), k = 1, . . . , k^ - 2, km , 

X^^'^km - I) = X^^'^km - I) + I. 

All the other coordinates remain unchanged, that is, 

x^'^ {k) = x^'^ {k) , j = l,...,d, k = l,...,ki, i = l,...,d, i ^ m . 

The sequence x is obtained from x by decreasing by one unit, the last jump of the 
coordinate x"^'^, that is x"^'^{km) — x™''^{km — 1). (Therefore the jump x"^'^{km — 1) — 
x"^'^{km — 2) is increased by one unit.) Denote by A'j.^x the number of good cyclical 
permutations of x. We claim that 

(4.20) > . 

To achieve this aim, first observe that k is a solution of the system (r, x), that is 

d 
i=l 

and note the straightforward inequality, 

x'''^ {k) > x^'^ {k) , i,j = l,...,d, k = l,...,ki. 

Therefore, if k is the smallest solution of (r, x) then it is also the smallest solution of (r, x). 
Moreover, for q = {qi, . . . ,qd) G N*^ satisfying < km — 1, one easily checks that the 
same inequality holds, that is 

•^qfk(^) — ^qfk(^) ' h j = ■ ■ ■ i d , k = 1, . . . , ki , 

SO that if k is the smallest solution of (r, Xq k) (i.e. Xq^k is a good cyclical permutation of 
x) then it is also the smallest solution of (r,Xq^k)- 

Now it remains to study the case where Qm = km — 1- Assume that Xq^k is a good 
cyclical permutation of x, but that Xq k is not a good cyclical permutation of x. Then in 
order to obtain the inequality (14.201) . we have to find a good cyclical permutation Xi^k of x 
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such that Zi^k is not a good cydical permutation of x. Let us define the sequence x G Sa 

"q,k ' 



with length k, which is obtained by decreasing by one unit the first coordinate of x^™''' 



that is 



where ei = (1, 0, ■ ■ ■ , 0), is the d dimensional unit vector. Since k is the smallest solution 
of (r, Xq^k), then k is the smallest solution of (r + ei,x) by the definition of x. Moreover, 
from Lemma 2.2 (ii), the system (ei,x) admits a smallest solution which is less than k. 
Let us call p this solution. Then p > and from Lemma 2.2 (ii), k — p is the smallest 
solution of (r, Xp^k)- 

Then let us consider the cyclical permutation of Xq k at p. It is a cyclical permutation 
of X that we shall denote by Xi^k- Note that x and Xq^k only differ from the last jump of 
x*"'^ and x^k^, more specifically, Xq^k(^m) — 5^q,k(^m ~ 1) = x{km) — x{km — 1) + 1. Then 
from the above constructions, we can see that Xi.k is obtained as follows: 

Xpk(/c), k < ki - Pi, i m, 



x^k) 



•'p,k 

^p,k(^i ~ Pi) + - {ki - Pi)), ki - Pi < k < ki, i m, 

•^Pjk (^"1 ~ Pm) + x'-™"-' {k — [km — Pm)) + Gil ~ Pm ^ k ^ k^ ■ 



Since k — p is the smallest solution of (r, Xp_k) and since Xi^k is strictly greater than Xp^k at 
point k — p, there is no solution to (r, Xi^k); on the (multidimensional) interval [0, k — p]. 
Moreover, also from the construction of (r, Xi^k); since p is the smallest solution of (ei, x), 
the only solution of (r, Xi^k) on the interval [k — p, k] is k. Therefore, the smallest solution 
of (r, Xi^k) is k. 

On the other hand, note the following identity 

(■^i,k(^)' h j — ■ ■ ■ k = 1, . . . , ki — Pi) = (ai'p^k(^)' h j — ^1 ■ ■ ■ k = 1, . . . , ki — Pi) , 

which can be seen directly from the definition of x. It follows that k — p is the smallest 
solution of the system (r, Xi_k)- But since p > 0, the sequence Xi_k is not a good permutation 
of X, and the inequality fl4.20p is proved. 

Let q = (0, . . . , 0, fcm — 1, . . . , 0), where km — ^ is the m-th coordinate of q and set 
y '■= 5:q,k, then by applying the same arguments as above to the chain y, we obtain that 

with obvious notations. But by reiterating km times this operation, we obtain again the 
chain x. This shows that equality holds in (I4.20p . that is A^^,* = Nj-^^- 

Finally, let ^ G 5*^ be a chain with length k, such that z^'^{ki) = x'^'^{ki), for all 
i,j = 1, . . . ,d, and such that for all i = 1, . . . , d, the sequences z^'\ j = 1, . . .d have at 
most one positive jump that occurs at the same point in {0, . . . , ki}. Then it is easy to see 
that the chain x can be obtained after several cyclical permutations and iterations of the 
transformation (14.191) applied to z, at any coordinate. Therefore A^j. ^ = A^r,x- Assume that 
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there is a good cyclical permutation z' of z. Then (r, z') is a simple system. Therefore, 
from Lemma [4.41 and Lemma [4. 5[ Ny-^^^t = N^^z = ^r,a; = det(a;*'-'(A;i)). □ 

Lemma 4.7. Let x G 5*^ with finite length and let x he the sequence of S^^^ which is associ- 
ated to X, as in fl2.8p . Then x and x have the same number of good cyclical permutations. 

Proof. Let k and n be the respective lengths of x and x. In particular, we have 
Hi = rj:^ , i = 1, . . . , d. Let q < k — 1^ be such that Xq^k is a good cyclical permutation of 

X. Then clearly, there is p = (pi, . . . ,pd) < n — 1^^, such that pi = Tqf . Set y = Xp^^ and 
let us check that 

(4.21) y = Xq,k . 

Define 9^'^\m) = min{u : x'^^'^^.{v) = —m} = minl^; : y^'^{v) = —m} and let rj^^^. be the 
cyclical permutation of the sequence r'-*-', as defined in fl3.13p . Then from the construction 
of y, we can check that 

O^Hm) = Tl%^{m) , m<h, 

from which we derive (I4.2ip . Moreover, since y is a good cyclical permutation of x, we 
deduce from Lemma 12.61 that y is a good cyclical permutation of x. 

Conversely, let q < n — 1,^ such that Xq^n is a good cyclical permutation of x. Then from 
part (z) of Lemma l2.2[ we must have rij = min{n : x*'*„.(n) = —ki}. Therefore, there 

exists p < k, such that qi = Tp} . Again, by setting y = Xq^n, we check that y = Xp^k and 
we deduce from Lemma [2.61 that y is a good cyclical permutation of x. □ 

Then we end the proof of Lemma 13.31 Let x G 5^ be as in this lemma, that is x has 
finite length n = (rii, . . . ,77,^) G N'^ and n is a solution of the system (r, x), where r = 
(ri, . . . , r^) G Z"^ is such that r > 0. Let A;- = — mino<n<ni 2;^* and set n- = t^K Then 

i 

n' < n, so that we can define the cyclical permutation y = Xn',n of x. By construction, 
the sequence y has length k, where ki = — x*'*(nj) and y'^'^{ki) = x*'-'(nj), i,j = 1, . . . ,d. 
Note that ki > 1, from the assumption x*'*(nj) 7^ 0. Moreover, k is a solution of the 
system (r,y). So, thanks to Lemma |4.6[ the number of good cyclical permutations of 
y is det(x*'-' (nj)) and from Lemmas 12.61 and 14.71 this is also the number of good cyclical 
permutations of y, the latter being clearly the number of good cyclical permutations 
of X. □ 
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